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A detector undergoing uniform acceleration a in a vacuum field responds just as 
though it were immersed in thermal radiation of temperature T = ha/2irkc. A 
simple, intuitive derivation of this result is given for the case of a scalar field in one 
spatial dimension. The approach is then extended to treat the case where the field 
seen by the accelerated observer is a spin- 1/2 Dirac field. 

I. INTRODUCTION 

Hawking 1 predicted that a black hole should radiate with a temperature T = hg/27rkc, 
where g is the gravitational acceleration at the surface of the black hole, k is Boltzmann's 
constant, and c is the speed of light. This results from the effect of the strong gravitation 
on the vacuum field. Shortly thereafter it was shown separately by Davies and Unruh that a 
uniformly accelerated detector in vacuum responds just as though it were in a thermal field 
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of temperature 2 ' 3 ' 4 ' 5 ' 6 



where a is the acceleration in the instantaneous rest frame of the detector. These results 
suggest profound consequences for the merger of quantum field theory and general relativity 
and sparked intense debates over unresolved questions that are still actively investigated 
today: (1) if black holes are not really "black," are naked singularities the ultimate fate of 
black holes, or will a long-sought fusion of quantum mechanics and general relativity into a 
coherent theory of quantum gravity prevent such occurrences?; (2) if a quantum mechanical 
pure state is dropped into a black hole and pure thermal (uncorrelated) radiation results, 
how does one explain the apparent non-unitary evolution of a pure state to a mixed state? 

This intriguing result of quantum field theory has arguably not been derived in any 
physically intuitive way. Numerous explicit and detailed calculations have appeared in the 
scientific literature over the last, roughly 30 years for a wide host of spacetimes. However, 
even for the "simplest" calculation involving a scalar field, the intricacies of field theory 
techniques, coupled with a forest of special function properties, makes most derivations 
intractable for the curious, casual scientific non-specialist. An investigation of the Unruh 
effect for the case of Dirac particles of spin 1/2 brings in a whole host of new machinery, least 
of which is the formulation of the Dirac equation in curvilinear coordinates (i.e. essentially in 
curved spacetime). This quickly goes beyond the expertise of most casual scientific readers. 
However, in both cases, this beautiful and important result can be stated quite simply: for 
a scalar field (bosons) the accelerated observer sees a Bose-Einstein (BE) distribution at 
temperature T given by Eq.(l), while for a spin-1/2 field the accelerated observer sees a 
Fermi-Dirac distribution at the same temperature. 

It is the purpose of this paper to present simplified derivations of Eq.(l) in a way that 
is suitable for advanced undergraduate or beginning graduate students and that elucidates 
the essential underlying physics of the Unruh effect. Once one accepts the simplest features 
of a quantized vacuum field, the result Eq.(l) emerges as a consequence of time dependent 
Doppler shifts in the field seen by the accelerated observer. 

In the following section the essential features of uniform acceleration for our purposes 
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are reviewed, and in Section 3 we use these results to obtain Eq.(l) in an almost trivial 
way based on the Doppler effect. In Section 4 this simple approach to the derivation of 
the temperature Eq.(l) is developed in more detail. In Section 5 we extend the previous 
calculations for scalar fields to spin-1/2 Dirac fields. We close with a brief summary and 
discussion. 



II. UNIFORM ACCELERATION 



Uniform acceleration is defined as a constant acceleration a in an instantaneous inertial 
frame in which the observer is at rest. The acceleration dv/dt in the lab frame is related to 
a by the Lorentz transformation formula 

,2\ 3/2 



a[l~- . (2) 



dv ( v 

Tt =a [ 1 -6, 

Integrating, and taking v = at t = 0, we have v(t) = at/yjl + a 2 t 2 /c 2 . The relation 
dt = dr I \j\ — v 2 / c 2 between the lab time (t) and the proper time (r) for the accelerated 
observer gives t(r) = (c/a) sinh(ar/c) if we take t(r = 0) = 0. The velocity v of the 
accelerated observer as seen from the lab frame can be expressed in terms of the proper time 

r as 

v (t) = ctanh (—) . (3) 



c / 

A straightforward integration of the above equations yields the well known hyperbolic orbit 
of the accelerated (Rindler) observer in the z direction 7 : 

t(r) = - sinh (—'] , z(r) = — cosh (—} ■ (4) 
a \ c J a \ c J 

Throughout this work we consider a > 0, i.e. the observer accelerates in the z direction. 



III. INDICATION OF THERMAL EFFECT OF ACCELERATION 

Consider now a plane-wave field of frequency ojk and wave vector K parallel or anti- 
parallel to the direction z along which the observer is accelerated. In the instantaneous rest 
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frame of the observer the frequency u' K of this field is given by the Lorentz transformation 
formula 

u k -Kv(t) c^[l -tanh(f )] 



= u K e- aT / c (K = +u K /c) (5) 



^\-v\t)I& ^l-tanh 2 (f) 

for K = +uk/c, i.e., for plane-wave propagation along the direction z of the observer's 
acceleration. For propagation in the — z direction, similarly, 



u' k (t) = uj K e aT ' c (K = -u K /c) . 



(6) 



Note that, for small values of ar, u' K = lvk(1 =F ot/c), the familiar Doppler shift. Equations 
(5) and (6) involve time- dependent Doppler shifts seen by the accelerated observer. 

Because of these Doppler shifts our accelerated observer sees waves with a time-dependent 
phase (p(r) = j T uj' K {T')dr' = (ukc/cl) exp(ar/c). We suppose therefore that, for a wave 
propagating in the — z direction, for which j T ijj' K {j')dr' = (cukc/cl) exp(ar/c), he sees a 
frequency spectrum S(Q) proportional to 



/oo , 
dre inT e i ( u} Kc/a)e aT > c 
-oo 



(7) 



Changing variables to y — e aT//c , we have 



d Te i^r e i(ui K c/a)e aT ' c 



a Jo 



dyy 



(iflc/a-1) i(uj K c/a)y 



c / iQc\ (u) K c 



where V is the gamma function 8 . Then, since 9 

/oo , 
dTe XlT e i(u, K c/a)e°T/ 
-oo 



a ) V a 

2 



-iQc/a 



(8) 



r (iOs) I = n/[(Qc/a) sinh(vrfic/a)], 
2 2ttc 1 



(9) 



tta e 27rnc / a - 1 ' 

The time-dependent Doppler shift seen by the accelerated observer therefore leads to the 
Planck factor ( e nn / kT — l)^ 1 with T = ha/2irkc, which is just equation Eq.(l). We obtain 
the same result in the case of a wave propagating in the +z direction. 

Note that the time-dependent phase can also be obtained directly by considering the 
standard nonaccelerated Minkowski plane wave exp[i</?±] = exp[i(Kz ± uxt)] and using 
equations Eq.(4): v 9 ±( r ) — Kz(t) ±u K t(r) = {uj K c/a) exp(±ar/c) with K = u K /c 10,11 . 
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IV. A MORE FORMAL DERIVATION 

The "derivation" of the temperature Eq.(l) just given leaves much to be desired. We 
have restricted ourselves to a single field frequency u>k, whereas a quantum field in vacuum 
has components at all frequencies. Moreover we have noted the appearance of the Planck 
factor but have not actually compared our result to that appropriate to an observer at rest 
in a thermal field. 

To rectify these deficiencies, let us consider a massless scalar field in one spatial dimension 
(z), quantized in a box of volume V 12 : 

0=1: ) l^ e ~ tujKt + • ( io ) 

Here K = ±uk/c, and clk and a) K are respectively the annihilation and creation operators 
for mode K ([a^,a^,] = Skk', [o>k,Q>k'] = 0). We use a caret ( * ) to denote quantum- 
mechanical operators. The expectation value ((dcj) / dt) 2 ) / 4nc 2 of the energy density of this 
field is V^ 1 J2k ^k[(o^k^k) + 1/2]- For simplicity we consider the field at a particular point 
in space (say, z — 0), since spatial variations of the field will be of no consequence for our 
purposes. 

For a thermal state the number operator a) K aK has the expectation value (e huJK / kT — l) -1 . 
Consider the Fourier transform operator 

1 r°° - f?irfir 2 \ 1/2 

The expectation value {g^(£l)g(£l')) in thermal equilibrium is therefore 



(12) 



We go to the limit where the volume of our quantization box becomes very large, V — > oo, so 
that we can replace in the usual fashion the sum over K by an integral: J2k ~^ (V/^) I dK 13 . 
Thus 

roo i I 

(tf(n)g(n')) = he 2 J ^ dK— ehn/kT — s(\k\c- n)S(n - of) 
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2hc/n 



-5(Q-n') . 



(13) 



e nn/kT _ i 

Now let us consider an observer in uniform acceleration in the quantized vacuum field. 
This observer sees each field frequency Doppler-shifted according to Eq.(5) and Eq.(6), and 
so for him the operator g(Q) has the form 



g(Q) = — / dTe inr Y 

yV ' 27T J-oo Y V "kV J 



(2Tthc 2 \ 1/2 



2lX J-oo x \ UJkV J 

j r (J K e~ i{ - eKU)KC l a)e ~ €KbJKT,c ] 



{e K uj K c/a)e-^K^/c 



a K e 



(14) 



where = \K\/K. Since a^|vacuum) = 0, only the a) K terms in this expression contribute 
to the vacuum expectation value (g^ {VL)g{Q!)) . Performing the integrals over r as before, 
and using {a K a) K ,) = 5kk>, we obtain 



2txo) \ V 



-ircQ/a 



E 



\ a 



ieif(f2— Q')c/a 



(15) 



K U K 

where we use the fact that the sum over k vanishes unless Q — Q'. In fact we show in the 

-2\ 



Appendix that the sum over K is (2Va/c )S(fl — f2'), so that 



na 



-5(0, - n') , 



which is identical to the thermal result Eq.(13) if we define the temperature by equation 
Eq.(l). 

Note that the expectation value {<xko} k ,) = Skk' involves the creation and annihilation 
operators of the accelerated observer and is taken with respect to the accelerated observer's 
vacuum, which is different from the vacuum seen by the nonaccelerated observer. This point 
is discussed more fully in Section VI. 



V. FERMI-DIRAC STATISTICS FOR DIRAC PARTICLES 



In the above we have considered a scalar field and have derived the Planck factor ( e m / kT — 
1) _1 indicative of Bose-Einstein (BE) statistics. We began with the standard plane- wave 
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solutions of the form exp[i(Kz ± to K t)] for the nonaccelerated Minkowski observer, and 
considered the time-dependent Doppler shifts as seen by the accelerated observer. For the 
spin-1/2 Dirac particles one would expect an analogous derivation to reproduce the Planck 
factor (e nn / kT + indicative of Fermi-Dirac (FD) statistics. 

We now show that this is indeed the case. Mathematically, the essential point involves the 
replacement iQc/a — > iVtc/a + 1/2 in the integrals in Eq.(7)- Eq.(9) 14 , and the relationship 
\T(iVLc/a + 1/2)| 2 = n/ cosh(7rf2c/a) 9 . Physically, this replacment arises from the additional 
spinor nature of the Dirac wave function over that of the scalar plane wave. In the case of a 
scalar field, only the phase had to be instantaneously Lorentz-transformed to the comoving 
frame of the accelerated observer. For non-zero spin, the spinor structure of the particles 
must also be transformed 15 , or " Fermi- Walker transported" 16 along a particle's trajectory 
to ensure that it does not "rotate" as it travels along the accelerated trajectory. This leads 
to a time-dependent Lorentz transformation of the Dirac bispinor of the form 17 S(t) = 
exp(7°7 3 ar/2c) = cosh(ar/2c) +7°7 3 sinh(ar/2c), where the 4x4 constant Dirac matrices 
are given by 



and a z = diagonal(l, —1) is the usual 2x2 Pauli spin matrix in the z direction. Acting on 
a spin up state ||) = [1,0, 1,0] T18 , S(t) gives S(t) ||) = exp(ar/2c) ||). Thus, for the spin 
up Dirac particle we should replace the plane wave scalar "wave function" exp[^(r)] used 
in Eq.(7) with exp(ar/2c) exp[i<^(r)] 19 . This leads to the replacement iQc/a — > iVtc/a + 1/2 
in Eq.(8), and therefore the result 



Comparing Eq.(17) with Eq.(9) we note the change of sign in the denominator from —1 for 
BE statistics to +1 for FD statistics. We also note that the prefactor in Eq.(9) involves the 
dimensionless frequency Vtc/a while in Eq.(17) the prefactor involves the factor uxc/a (the 
argument of the exponential in the distribution function is still TiVt/kT with the same Unruh 
temperature T = ha/27ikc). This is no cause for concern, since in fact a single Minkowski 
frequency ujk is actually spread over a continuous range of accelerated (Rindler) frequencies 





(17) 
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Q with peak centered at Q = w K 20 . This allows us to replace uo K by Q in the final result. 
(This is explicitly evidenced by the delta function 5(oo K — Q) in Eq.(ll) - Eq.(16) in the 
comparison of the thermal and accelerated correlation functions.) 

For the case of the spin up Dirac particle, the more formal field-theoretic derivation of 
Section IV proceeds in exactly the same fashion, with the modification of the accelerated 
wave function from exp[i<^(r)] — ► exp(ar/2c) exp[i<^(r)] and the use of anti-commutators 
{aKyC^K'} = $kk' for the quantum-mechanical creation and annihilation operators instead 
of the commutators appropriate for scalar BE particles. For the correlation function we find 

= ^f-.m - no , (is) 

the FD analogue of Eq.(16). 



VI. SUMMARY AND DISCUSSION 

In the usual derivation of the Unruh temperature Eq.(l), 2 ' 3 ' 4 ' 5 ' 6 one solves the wave (or 
Dirac) equation for the field mode functions in the Rindler coordinates Eq.(4) and then 
quantizes them. Since the hyperbolic orbit of the accelerated observer is confined to the 
(right Rindler) wedge z > \t\ (with "mirror" orbits in the left Rindler wedge z < —\t\) 
it turns out that the vacuum seen by the accelerated observer in say, the right Rindler 
wedge is different than the usual Minkowski vacuum (defined for all z and t) seen by the 
unaccelerated observer. The inequivalence of these vacua (and hence the Minkowski vs. 
Rindler quantization procedures 21 ) is due to the fact that the right and left Rindler wedges 
are causally disconnected from each other. This can be easily seen by drawing a Minkowski 
diagram in (z, t) coordinates and observing that light rays at ±45° emanating from one 
wedge do not penetrate the other wedge. Hence the Minkowski vacuum that the accelerated 
observer moves through appears to him as an excited state containing particles, and not as 
the vacuum appropriate for the right Rindler wedge. The Bose-Einstein distribution with 
Unruh temperature T for scalar fields (Fermi-Dirac for Dirac fields) then arises by considering 
the expectation value of the number operator a} R aR for the accelerated observer (in the 
right Rindler wedge) in the unaccelerated Minkowski vacuum \0m), i-e. (0m\ 
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(exp(hQ I kT) ± (with upper sign for scalar fields and lower sign for Dirac fields). This 
is referred to as the Thermalization Theorem by Takagi 5 . 

In this work we take a slightly different viewpoint 22 . For a scalar field, we first consider an 
unaccelerated Minkowski observer in a thermal state and find that the expectation value of a 
field correlation function is proportional to the Bose-Einstein distribution. We then consider 
the calculation of this correlation function again, but this time for an accelerated Rindler 
observer in his Rindler vacuum 23 state |0^) , such that for a single mode, (0^|ci^ajj|0^) = 1. 
The new feature now is that from his local stationary perspective, the accelerated observer 
perceives all Minkowski frequencies (arising from the the usual plane waves associated with 
Minkowski states) as time-dependent Doppler shifted frequencies. 

The derivation presented in this work shows why quantum field fluctuations in the vacuum 
state are crucial for the thermal effect of acceleration: (g^(Q)g(Q')) is nonvanishing because 
the vacuum expectation (a^-a^) 7^ 0. But there's more to it than that, because (oR-a^) is 
also nonvanishing for an observer with a = 0. For such an observer, however, 



for the case of scalar particles. In other words, the thermal effect of acceleration in our 
model arises because of the nontrivial nature of the quantum vacuum and the time-dependent 
Doppler shifts seen by the accelerated observer. For the case of Dirac particles, the essential 
new feature is the additional spinor structure of the wave function over that of the scalar 
plane wave. In order to keep the spin "non-rotating" in the comoving frame of the acceler- 
ated observer, the Dirac bispinor must be Fermi- Walker transported along the accelerated 
trajectory, resulting in an additional time- dependent Lorentz transformation. Formally, this 
induces a shifting of iQc/a — > iQc/a + 1/2 in the calculation of relevant gamma function-like 
integrals, leading to the FD Planck factor. 

In the following we briefly discuss the relationship of our correlation function to those used 
in the standard literature on this subject and point out a not widely appreciated curious 
subtlety relating details of the spatial Rindler mode functions (which we have ignored in our 
model) to the statistics of the noise spectrum seen by the accelerated observer. 

In our model, we have not given any motivation for using the correlation function 




(19) 
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(g^(Q)g(Q')) aside from the fact that we could calculate it for a nonaccelerated observer 
in a thermal field and for a uniformly accelerated observer in vacuum and compare the 
results. It is easy to show that a harmonic oscillator with frequency uj and dissipation 
coefficient 7, linearly coupled to the field Eq.(10), reaches a steady-state energy expectation 
value 

{E) K r dn r gg> . . , (20) 

Jo Jo (fi — ivo — -ou + 17) 

which offers some motivation for considering (gt(f2)g(f2')). In fact, it can be shown that 
(E) = [e hu} °/ kT - l]" 1 , which shows again that our accelerated observer acquires the charac- 
teristics appropriate to his being in a thermal field at the temperature T = ha/2nkc. 

In an extensive review of the Unruh effect, Takagi 5 (Chapter 4) utilizes the quantum 
two-point correlation (Wightman) function gw(j,T r ) = (0m | 4>{t) <ft {t') \ 0m) to determine 
the power spectrum of the vacuum noise seen by the accelerated observer for a scalar field 
via 

5(n)=lim/ e-^-^W(r), 

which is very much in the spirit of our calculation. Here the field operator 0(r) is ex- 
panded in terms of the Rindler mode functions and involve creation and annihilation op- 
erators for both the right and left Rindler wedges. Takagi shows the remarkable, though 
not widely known result, that for a scalar field in a Rindler spacetime of dimension n, 
S n (£l) ~ f n (Q)/[exp(hQ/kT) — (—1)™]. For even- dimensional spacetimes (e.g. n = 2 consid- 
ered in this work, or the usual n — 4) S n (Q) is proportional to the Bose-Einstein distribution 
function [exp(Ml / "kT) — l] -1 , and essentially reproduces our equation Eq.(16) (up to powers 
of flc/a, contained in the function f n (Q)). However, for n odd, S n (Q) is actually propor- 
tional to the Fermi-Dirac distribution [exp(hQ/kT) + l] -1 . For the case of Dirac particles 
the opposite is true, namely for even spacetime dimensions S n (Q) is proportional to the FD 
distribution and for odd spacetime dimensions S n (Q) is proportional to the BE distribution. 
This curious fact arises from the dependence of S n (£l) on two factors in its calculation. The 
first is the above mentioned Thermalization Theorem, namely that the number spectrum of 
accelerated (Rindler) particles in the usual nonaccelerated Minkowski vacuum is proportional 
to the BE distribution function for scalar fields and is proportional to the FD distribution 
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function for Dirac fields. The second factor that switches the form of S n (Q) from BE to 
FD depends on the detailed form of the Rindler mode functions (see Takagi 5 , Chapter 4 for 
more complete details). Though the trajectory of the accelerated observer takes place in 
"1+1" dimensions (i.e. the (z,t) plane), the quantum field exists in the full n-dimensional 
spacetime, and thus S n (Q) ultimately depends of the form of the mode functions in the full 
spacetime. In spacetimes of even dimensions the number spectrum of Rindler particles in 
the Minkowski vacuum and the noise spectrum of the vacuum fluctuations (i.e., the response 
of the accelerated "particle detector") both depend on the same distribution function, and 
these two effects are often incorrectly equated. 

In our simplified derivation we have bypassed this technicality by performing our calcu- 
lations in "1+1" dimensions (i.e., n = 2). We have concentrated on the power spectrum of 
vacuum fluctuations as seen by a particle detector carried by the accelerated observer. We 
have shown that in "1+1" dimensions the spectrum of fluctuations is proportional to the 
Bose-Einstein distribution function for scalar fields and to the Fermi-Dirac distribution for 
spin-1/2 fields, with the Hawking-Unruh temperature defined by Eq.(l). The dependence 
of the noise spectrum on these distribution functions is ultimately traced back to the time- 
dependent Doppler shifts as seen by the accelerated observer as he moves through the usual 
nonaccelerated Minkowski vacuum. It is hoped that the simple calculations exhibited here 
are straightforward enough to give an intuitive understanding of the essential physical origin 
of the Hawking-Unruh temperature experienced by a uniformly accelerated observer. 

Exercise for the student Discuss when the Hawking-Unruh temperature from Eq.(l) 
would become physically detectable by utilizing the expression a = GM/r 2 for the 
gravitational acceleration of a test mass at a distance r from a mass M, and determining T 
at the surface of the earth, the sun and a Schwarzschild black hole. 



Appendix 

Converting the sum over K to an integral, we have, for K > 0, 

y, 1 fuj K C\ ieK ^-^ c / a _ V f°° dK 1 fUKC\ i{n ~ n ' )c,a . . 

u K \ a J 2tt Jo u K \ a J 
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Letting x = \og{ujKc/ a), we can write this as 



2-kc 



I" dxe- lx{n - n ' )c/a = ^5(n - Q') . (22) 

J-oo C 2 



The same result is obtained for the sum over K < 0, so that the sum over all K is 
{2Va/c 2 )5{Q-Q'). 
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Phys. Rev. D22, 1935 (1980). 

19 The spinor Lorentz transformation S(t) does not mix spin components. Thus, for e.g., a spin 
up Minkowski state remains a spin up accelerated (Rindler) state. We can therefore drop the 
constant spinor ||) from our calculations and retain the essential, new time-dependent 
modification exp(ar/2c) to the plane wave for our Dirac "wave function." 

20 Reference [5], Section 2, in particular Eqs. (2.7.4) and (2.8.8). 

21 S.A. Fulling, Phys. Rev. D 7, 2850 (1973); Aspects of Quantum Field Theory in Curved 
Space-Time, Cambridge Univ. Press, N.Y. (1989). 

22 This viewpoint is also taken in a different derivation of the Unruh-Davies effect in PW. 
Milonni, The Quantum Vacuum, (Academic Press, NY, 1994), p60-64. 

23 The unaccelerated Minkowski vacuum \0m) is unitarily related to the Rindler vacuum 
|0_r) <8> \0l) via \0m) = S(r)\0n} <8> \0l) where S(r) is the squeezing operator 

S(r) = exp[r(di?dL — a^ci^)]. Here the subscripts R and L denote the right (z > \t\) and left 
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(z < —\t\) Rindler wedges, respectively, and |0#) is the Fock state of zero particles in the right 
Rindler wedge and |0l) is the Fock state of zero particles in the left Rindler wedge. Since these 
right and left Rindler wedges of Minkowski spacetime are causally disconnected from 
eachother, the creation and annihilation operators a R ,a,R and a^,d^ live in the right and left 
wedges respectively and mutually commute with with each other, i.e. [or, ajj = 0, etc. . . . 
Since, physical states that live in the right wedge have zero support in the left wedge (and visa 
versa), they are described by functions solely of the operators a R ,a,R appropriate for the right 
wedge, i.e. \^ r) = f(an,a R ) |0#) <8> |0^) = \i/)r) <8> It is i n this sense that we can speak of 
| Or) as the "vacuum" for the right Rindler wedge, and similarly |0l) as the "vacuum" for the 
left Rindler wedge. 



